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ABSTRACT 

OfFshore  jets  or  filaments  are  a  well  observed  phenomenon  in  the  California 
Current  System.  Their  surface  and  subsurface  structures  as  well  as  their  spatial  and 
temporal  variability  are  well  documented.  The  methods  by  which  these  filaments  are 
formed  is  not  as  well  understood  as  their  structural  form.  The  most  common 
generation  theories  are:  (1)  variations  in  wind  stress  coupled  with  topographic 
irregularities,  (2)  dynamic  instability  and  (3)  geostrophic  turbulence.  In  this  study  an 
attempt  is  made  to  identify  the  factor(s)  responsible  for  filament  formation  south  of 
Cape  Mendocino  by  numerically  simulating  the  California  Current  System  using  a  two- 
layer,  non-linear,  primitive  equation  model.  It  is  shown  that  baroclinic  instability  is 
the  primary  method  by  which  filaments  are  formed  in  this  region  and  that  other  factors 
such  as  barotropic  instability,  friction  and  non-linearities  can  alter  the  characteristics  of 
the  fastest  growing  baroclinically  unstable  wave. 
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I.  INTRODUCTION 

A,       THE  CALIFORNIA  CURRENT  SYSTEM 

The  California  Current  System,  CCS,  has  been  extensively  studied  by  Ilickey 
(1979).  The  CCS  consists  of  four  separate  currents  that  are  spatially  and  seasonally 
distributed  off  the  West  Coast  of  the  United  States.  They  are:  (1)  the  California 
Current,  a  southward  flowing  surface  jet;  (2)  the  California  Undercurrent,  a  northward 
flowing  sub-surface  jet;  (3)  the  Davidson  Current,  a  poleward  flowing  surface  jet  north 
of  Point  Conception;  and  (4)  the  Southern  California  Countcrcurrent,  a  northward 
flowing  surface  current  south  of  Point  Conception  in  the  California  Bight  (Hickcy, 
1979).  Figure  1.1  shows  the  average  seasonal  configurations  from  Ikeda  and  F-mcry 
(1984)  for  the  CCS.  Typical  surface  jet  velocities  arc  on  the  order  of  20  to  40  cm/s. 
The  California  Undercurrent  has  a  typical  strength  of  about  5-10  cm  s.  For  the 
profiles  on  the  left,  where  the  California  Undercurrent  is  present,  the  surface  jet  depth 
is  about  150  m  while  for  the  profiles  without  a  subsurface  jet,  the  surface  jet  extends 
down  to  about  300  m. 
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Figure  1.1     Seasonal  Configurations  of  the  California  Current 
System  (Ikeda  and  Emery,  1984). 


B.       FILAMENTS  IN  THE  CALIFORNIA  CURRENT  SYSTEM 

The   presence   of  offshore  jets   or  filaments   in  the   CCS   is   well   documented. 
Satellite  radiometry  and  color  scanning  along  with  intensive  in  situ  data  such  as  that 


obtained  during  the  Coastal  Ocean  Dynamics  Experiment  (CODE)  have  provided  large 
data  bases  for  these  features. 

From  an  observational  standpoint,  filaments  are  well  documented.  Definitive 
studies  by  Flament  et  al.  (1985),  Kosro  (1987)  and  Kosro  and  Huyer  (1986),  for 
example,  have  greatly  expanded  our  knowledge  of  filament  surface  and  subsurface 
structures  along  with  their  spatial  and  temporal  variability. 

Filaments  consist  of  strong  baroclinic  jets  flowing  offshore  at  peak  velocities  of 
0.5  m  s  or  greater  (Rienecker  ei  al.,  1985),  with  narrow  return  fiows  to  the  south 
(Flament  e[  al.,  1985).  The  offshore  extent  of  filaments  is  typically  on  the  order  of  200 
to  300  km  \\ith  maxima  of  around  500  km  (Brink  and  Hartwig,  1985).  Filaments 
advect  cold  upwelled  coastal  waters  away  from  the  continental  shelf 

Filament  widths  are  on  the  order  of  50  km  (Flament  et  al.,  1985).  The 
wavelength  associated  with  their  occurrence  at  numerous  points  along  the  coast  is 
roughly  100  to  400  km  for  the  region  to  the  north  of  Cape  Mendocino  (Ikeda  and 
Emer>'.  1984).  Formation  time  scales  are  generally  on  the  order  of  weeks  with  lifetimes 
on  the  order  of  a  month.  Figure  1.2  shows  the  surface  velocity  field  derived  from 
satellite  feature  tracking  for  a  filament  studied  by  Flament  et  al.  (1985). 

Filaments  are  seasonal  in  nature.  Few  have  been  observed  during  late  winter 
through  early  spring.  Their  first  appearance  is  off  the  Baja  Peninsula  in  March. 
Between  April  and  May  they  occur  between  Point  Conception  and  Cape  Mendocino 
and  by  late  summer  they  reach  their  maximum  northward  extent  off  the  Oregon  coast. 
Generally,  filaments  have  been  more  observed  south  of  Cape  Mendocino.  (Brink  and 
Hartwig,  1985) 

C.       FILAMENT  FORMATION 

The  precise  formation  process  for  filaments  is  not  as  well  understood  as  their 
characteristics  once  formed.  Numerous  theories  abound  with  the  most  prevalent 
causative  forces  being:  (1)  coastally  linked  processes  associated  with  variations  in  wind 
stress  coupled  with  irregularities  in  the  alongshore  topography,  (2)  dynamical 
instabilities  of  the  CCS,  and  (3)  geostrophic  mesoscale  eddy  field  interactions  with  the 
CCS.  It  is  possible  and  indeed  probable  that  more  than  one  theory  may  be  responsible 
for  filament  formation. 

Ikeda  and  Emery  (1984)  demonstrated  by  a  combination  of  observational  data 
and  numerical  experiments  that  coastal  irregularities  off  the  coast  of  Northern 
California  and  Oregon  are  strongly  correlated  with  filament  formation.    They  also 
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Figure  1.2    Surface  Velocity  Field  of  a  Filament  (Flament  ei  ai,  1985). 

showed  that  baroclinic  instability  was  probably  the  primary  causative  force  north  of 
Cape  Mendocino.  Batteen  (1987)  demonstrated  that  baroclinic  instability  could  occur 
in  the  CCS  due  to  imposing  an  unstable  jet  profile  similar  to  the  early  summer  profile 
shown  in  Figure  1.1.  Chen  (1986)  showed  that  baroclinic  instability  coupled  with 
topographic  variations  are  important  to  the  destabilization  of  the  CCS  in  the  vicinity 
of  Cape  Mendocino.  Mysak.  (1977)  also  demonstrated  the  importance  of  baroclinic 
instability  in  a  numerical  study  of  the  CCS  off  Vancouver  Island.  Mooers  and 
Robinson  (1984)  documented  an  eddy  pair  off  Northern  California,  while  Owen  (1980) 
and  Huyer  et  al.  (1984)  observed  strong  eddy-like  features  in  the  northeastern  Pacific 
close  to  the  continental  slope. 

D.       PURPOSE  OF  THIS  STUDY 

The  objective  of  this  study  is  to  determine  the  causative  factor(s)  involved  with 
filament  formation  along  the  coast  of  California  south  of  Cape  Mendocino.  A  two- 
layer,  non-linear,  primitive  equation  model  with  topography  is  used  to  simulate  the 
dynamics  of  the  CCS.  The  results  strongly  point  to  baroclinic  instability  as  the 
primary  generation  mechanism  for  filaments  in  this  region. 
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The  specifics  of  the  model  and  the  various  numerical  experiments  are  covered  in 
Chapter  II.  The  linear  baroclinic  instability  problem  is  formulated  in  Chapter  III  and 
the  results  applied  to  the  data  from  the  numerical  experiments  to  verify  the  claim  of 
baroclinic  instability  to  the  CCS.  In  Chapter  IV,  the  results  are  discussed  and 
conclusions  presented. 
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II.  NUMERICAL  SIMULATION  OF  THE  CALIFORNIA  CURRENT 

A.       THE  NUMERICAL  MODEL 

L  Model  Equations 

Simulation  of  the  CCS  is  performed  using  a  two-layer,  primitive  equation, 
semi-implicit  numerical  scheme  initially  developed  by  Hurlburt  (1974).  Ihis  scheme 
has  been  employed  in  numerous  ocean  mesoscalc  circulation  studies  (e.g.,  Hurlburt  and 
Thompson,  1980,  1982;  Smith  and  O'Brien,  1983;  and  Smith,  1986).  Linear  test  cases 
have  been  run  for  comparison  with  linear  analytic  solutions  to  show  that  this  approach 
is  valid  (Smith  and  Reid,  1982).  Motion  in  each  layer  is  governed  by  a  momentum 
equation  (2.1)  and  a  continuity  equation  (2.2). 


tit 


(2.1) 


— 1+   y.V.    =  0 


(2.2) 


Where   i  =  layer  index  (i=  1  upper,  i=  2  lower) 


Variables  and  notation  are  defined  in  the  Appendix.    The  iluid  is   assumed  to   be 
hydrostatic  and  Boussinesq.     The  density  in  each  layer  is  constant.     1  he  clfccts   o[ 
winds,  tides,  thermodynamics  and  thermohaline  mixing  are  not  included. 
2.  Model  Domain 

A  rectangular  region  (1100  km  x  800  km)  was  divided  into  10  km  by  10  km 
squares  to  form  the  grid  for  the  numerical  model  finite  diflerencing.  By  rotating  the 
grid  through  a  desired  angle,  the  model  can  be  used  to  simulate  the  ilow  within  the 
CCS  as  indicated  by  Figure  2.1. 

Topography  is  included  by  approximating  the  major  features  apparent  in  the 
Mendocino  Escarpment  region.  Figure  2.2  shows  the  digitized  topography  in  the 
vicinity  of  the  Mendocino  Escarpment  (Chen,  I98C)  and  Figure  2.3  shows  the 
topography  used  in  the  model.  Principal  features  of  the  topography  used  in  the  model 
are:    (1)  a  shelf  region  of  500  m  depth  extending  out  50  km  from  the  coastline  which  is 
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Figure  2.1     Geographical  Positioning  of  iModel  Domain. 
Solid  box  represents  a  basin  rotation  (y)  of  90"*,  dashed  box  for  rotation  of  120' 
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a  uniform  vertical  wall;  (2)  a  shelf  slope  region  extending  from  50  km  to  ~  100  km  with 
an  exponentially  decreasing  slope;  (3)  a  modeled  Mendocino  Escarpment  (without  a 
ridge)  which  separates  a  3000  m  depth  region  to  the  north  from  4500  m  water  depth  to 
the  south.   The  topography  rotates  with  the  model  domain. 
3.  Boundary  Conditions 

Using  the  orientation  of  Figure  2.1,  no-slip  boundary  conditions  are  set  on  the 
eastern  and  western  boundaries.  The  northern  boundary  condition  consists  of  a 
prescribed  inflow  in  the  upper  layer.  For  simulated  lower  layer  flows  to  the  south  or 
north  the  lower  layer  boundary  consists  of  either  a  prescribed  inflow  or  a  radiation 
condition  (Camerlengo  and  O'Brien,  1980)  respectively.  The  southern  boundary 
condition  is  set  up  just  opposite  the  northern.  The  upper  layer  always  has  a  radiating 
condition  and  the  lower  is  radiating  or  prescribed  for  southward  and  northward  lower 
layer  flows  respectively. 

Obviously,  these  boundary  conditions  do  not  reflect  totally  the  situation  in  the 
CCS.  By  putting  a  closed  boundary  to  the  west,  there  will  inevitably  be  shear 
produced  vorticity  that  is  not  reaUstic.  However,  it  is  believed  that  the  east-west 
dimension  of  the  domain  is  large  enough  so  that  the  vorticity  production  at  the  western 
boundar\'  will  not  affect  the  dynamics  of  the  CCS  which  is  located  near  the  eastern 
boundarv'. 

B.       NUMERICAL  EXPERIMENTS 

1.  The  Reference  State 

In  order  to  filter  through  the  many  possible  generating  mechanisms  for 
filaments,  variation  of  selected  parameters  is  necessary.  The  problem  is  twofold:  (1) 
selective  variation  of  parameters  may  cause  model  deviations  from  the  CCS 
configuration  such  that  extrapolation  of  results  from  the  model  to  the  CCS  is 
impossible.  The  results  would  still  be  important  from  a  dynamical  standpoint  but 
without  the  connection  with  the  'real  world'  their  usefulness  is  limited;  (2)  isolation  of 
specific  model  responses  by  changes  in  model  parameter(s)  can  be  inaccurate  due  to 
complex  cross  correlations  between  varied  and  non-varied  parameters. 

In  order  to  assess  the  model  output  variations  resulting  from  parameter 
changes,  a  reference  state  is  estabUshed.  Instead  of  a  simple  configuration  such  as  flat 
bottom  barotropic  flow,  the  reference  state  is  selected  to  be  as  complex  as  possible 
involving  such  things  as  baroclinicity,  topography  and  the  presence  of  a  pre-existing 
eddy  (all  suggested  to  be  dynamically  important  in  Chapter  I).    As  mentioned  before, 
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Figure  2.2    Actual  Topography  in  Vicinity  of  Mendocino  Escarpment  (Chen,  1980). 
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Figure  2.3     Smoothed  Topography  used  in  Numerical  Simulations. 
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applicability  of  the  results  to  the  CCS  is  vital  and  iise  of  a  reference  state  as  near  as 
possible  to  the  CCS  configuration  keeps  the  study  focus  on  CCS  dynamics  rather  than 
general  dynamics. 

The  reference  state  consists  of  a  two-layer  geostrophically  balanced,  Gaussian 

jet  of  the  form, 

h^(y)  =  A^[l-exp{-y2;2L2)]       i=  1.2  (2.3) 

where  L  is  the  e-folding  width  scale  of  the  jet  and  the  y-axis  oriented  east-west  as 
shown  in  Figure  2.1.  The  flow  strength  is  determined  by  the  amplitude  coefficient  A- 
which  determines  the  magnitude  of  the  surface  and  interface  height  anomalies  through 
hj(y).   Additionally,  a  two-layer  geostrophically  balanced,  axisymmetric,  Gaussian  eddy 

is  included, 

h-(x,y)  =  A^[l-exp(-R2/2L2)]    i=  1,2  (2.4) 

where  R^  =  (x^  +  y^)  and  represents  the  radial  distribution  of  the  eddy. 

Mean  upper  layer  thickness  is  set  at  150  m  (H^)  which  corresponds  to  a  first 
internal  Rossby  radius,  R^,  of  ~  18  km.  This  is  comparable  to  the  ~25  km  R^^  for  the 
CCS  (Brink  and  Hartwig,  1985).    H2  is  dependent  on  the  topography  depth,  D(x,y),  by, 

H2(x,y)  =  4500  -  Hj  -  D(x,y)  (2.5) 

The  jet  e-folding  scale,  L,  is  set  at  25  km.  The  jet  is  fixed  parallel  to  the  coast  over  the 
slope  region  and  the  eddy  positioned  offshore  from  the  jet.  The  spatial  relationships 
between  the  jet,  eddy  and  topography  are  indicated  in  Figure  2.4.  In  Figure  2.4  and  in 
all  subsequent  figures  of  model  output  fields,  tick  marks  along  the  x  and  y  axis 
represent  a  distance  of  275  km  and  200  km  respectively.  The  lines  o[  constant  surface 
height  anomaly  can  be  thought  of  as  streamlines  representing  the  barotropic  part  of 
the  current.  The  interface  anomaly  is  sloped  opposite  the  surface  with  an  amplitude 
such  that  its  baroclinic  component  in  the  opposite  direction  sets  up  the  selected  flow 
magnitude  in  the  lower  layer. 

The  initial  jet  and  eddy  strength  is  set  at  0.4  m/s  (southward  for  jet  and 
anticyclonic  for  the  eddy)  in  the  upper  layer  and  0.01  m/s  in  the  same  direction  in  the 
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Figure  2.4     Reference  State  Flows  and  Topography. 

Reference  state  initial  conditions  of  (a)  surface  hei2ht  anonialy(cm)  and  (b)  lieicht  of 
bottom  above  4500  m  (m).  Contour  intervals  arc"  (a)  2.5  cm  and  (b)  250  m.  lick 
marks  represent  275  km  m  the  x  direction  and  200  km  in  the  y  direction. 
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lower  layer.  The  horizontal  Laplacian  frictional  coefficient,  A^,  is  set  at  500  m^/s 
which  is  roughly  in  the  middle  of  the  range  of  A^  (100-1000  n?!s)  typically  used  in 
numerical  studies.  A  summary  of  model  parameters  for  the  reference  state  is  shown  in 
Table  1. 


TABLE  1 

■ 

REFERENCE  STATE  MODEL  PARAMETERS 

PAR.AMETER 

SYMBOL 

VALUE 

E-W  Basin  Extent 

Lv 

800  km 

N-S  Basin  Extent 

Lx 

1100  km 

Initial  Upper  Layer  Thickness 

Hi 

150  m 

Initial  Lower  Layer  Thickness 

H2(x,y) 

Variable 

Maximum  Basin  Depth              (Hj  4-  yi2^msLX 

4500  m 

Depth  of  Bottom 

D(x,y) 

Variable 

Coriolis  Parameter  (@  40*  N  lat 

)      ^0 

0.94  X   lO-'^s'^ 

df'dy 

Po 

2.0  X   10-11  j^-lg-l 

Gravitational  Acceleration 

g 

9.8  m  s-2 

Reduced  Gravity 

g' 

2.0  X   10-2  ^  5-2 

Horizontal  Eddy  Viscosity  Coef 

^h 

5.0  X   lO^m^s-1 

Time  Step 

At 

3600  s 

Grid  Increment  in  x-direction 

Ax 

10  km 

Grid  Increment  in  y-direction 

Ay 

10  km 

Jet  Maximum,  Upper  Layer 

^'ij 

0.4  m  s"    (southward) 

Jet  Maximum,  Lower  Layer 

^2j 

0.01  m  s-1  (southward) 

Eddy  Maximum,  Upper  Layer 

Vie 

0.4  m  s-1  (anticyclonic) 

Eddy  Maximum,  Lower  Layer 

Vie 

0.4  m  s'l  (anticyclonic) 

Jet  Position  Offshore 

H 

50  km  (from  east  wall) 

Eddy  Position  Offshore 

Le 

200  km  (from  east  wall) 

Basin  Rotation  Angle 

Y 

90"  (from  east) 

Figures  2.5  and  2.6  show  the  reference  state  model  output  fields  for  height 
anomalies  and  upper  layer  velocities  respectively.    For  this  and  all  subsequent  model 
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output  field  plots,  surface  height  anomaly  is  contoured  in  cm,  interface  height  anomaly 
is  contoured  in  m  and  the  velocities  are  contoured  in  cm,  s.  Based  on  the  directional 
characteristics  and  the  magnitudes  of  the  flows,  it  is  apparent  that  some  phenomenon 
very  similar  to  filament  formation  is  occurring.  Of  particular  note  is  the  velocity 
structure  of  the  offshore  flow  versus  the  inshore  flow.  The  offshore  flow  has  strong 
unidirectional  flow  while  the  inshore  flow  is  slower  and  has  more  directional  variability. 
The  maximum  offshore  velocity  by  day  30  is  about  50  cm/s  and  the  inshore  flow 
maximum  about  40  cm/s.  This  type  of  structure  and  magnitude  of  flow  correlates  well 
with  observed  data  in  filament  studies  (Flament  et  ai,  1985;  Kosro,  1987). 

There  are  also  two  important  signatures  of  baroclinic  instability  in  the  output 
fields:    (1)  wave  growth  in  the  initially  nearly  quiescent  lower  layer  is  indicative  of 
vertical  energy  transfer  from  the  upper  to  lower  layer,  and  (2)  there  is  a  180°  phase 
shift  between  the  waves  in  the  upper  layer  and  the  waves  in  the  lower  layer. 
2.  Variation  of  Parameters 

Structures  shown  in  Figure  2.5  are  consistent  with  observations  of  the  form 
and  magnitude  of  filaments.  In  order  to  isolate  potential  causes  for  filaments,  the 
assumptions  that  the  wave  structures  in  Figure  2.5  represent  filament  formation  and 
growth  are  made.  Since  the  reference  state  is  as  complex  as  possible,  the  variations  of 
the  numerical  simulations  are  largely  one  of  simplification  in  specific  areas  of  the  model 
to  see  whether  filaments  are  altered  or  made  non-existent. 

All  of  the  numerical  simulations  are  for  30  days  and  certain  parameters  from 
Table  1  are  never  varied.  Table  2  provides  a  Ust  of  those  parameters  which  are  varied 
in  the  different  experiments  and  what  their  reference  state  value  is. 

a.  Experiment  No.  1  {H,  Increased  to  250  m) 

Upper  layer  thickness  was  increased  to  250  m  to  determine  the  model 
sensitivity  to  the  vertical  structure.  Figure  2.7  shows  the  model  output  fields  for 
surface  and  interface  height  anomalies.  From  the  fields,  it  is  apparent  that  the  wave 
structures  are  significantly  different  than  those  in  the  reference  state.  They  are  weaker 
and  seem  to  stop  growing  halfway  through  the  sequence.  This  points  to  the 
importance  of  the  vertical  structure  to  the  wave  development  process. 

b.  Experiment  No.  2  {Eddy  Removed) 

The  eddy  was  removed  in  this  experiment  to  determine  whether  it  was 
responsible  for  the  wave  structures  by  its  interaction  with  the  coastal  jet.  Figure  2.8 
shows  the  output  fields  for  surface  and  interface  height  anomalies.    The  results  are  very 
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Figure  2.5     Model  Output  Fields  of  Height  Anomalies 
for  Reference  State. 

Height  anomaly  fields  for  (a)  surface  at  dav  15  (cm),  (b)  surface  at  day  30  (cm),  (c) 
interface  at  day  15  (m),  and  (d)  interface  at  day  30  (m). 
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(a) 


(b) 


Figure  2.6     Model  Output  Fields  of  Upper  Layer  Velocity 
for  Reference  State. 

Velocity  fields  (cm/s)  for  (a)  upper  layer  at  day  15  and  (b)  upper  layer  at  day  30. 
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TABLE  2 

VARIATION  OF  MODEL  PARAMETERS 

PARAMETER  VALUES  (Notation 

from  Table  1) 

Hi 

D(x,y) 

vii 

^2i 

Vie 

V2e 

H 

Le 

Y 

Exp.  p- 

ref. 

150 

-.4 

-.01 

-.4 

-.01 

50 

200 

90 

1. 

250 

-.4 

-.01 

-.4 

-.01 

50 

200 

90 

2. 

150 

-.4 

-.01 

ne 

ne 

50 

200 

90 

3. 

150 

-.4 

-.01 

-.4 

-.01 

150 

300 

90 

4. 

150 

0 

-.4 

-.01 

-.4 

-.01 

50 

200 

90 

5. 

150 

-.4 

-.05 

-.4 

-.05 

50 

200 

90 

6. 

150 

-.4 

.01 

-.4 

.01 

50 

200 

90 

7. 

150 

-.4 

-.01 

-.4 

-.01 

50 

200 

0 

8. 

150 

-.4 

-.01 

-.4 

-.01 

50 

200 

180 

9. 

150 

0 

-.4 

-.01 

ne 

ne 

50 

200 

90 

10. 

150 

-.4 

-.01 

ne 

ne 

50 

200 

120 

11. 

150 

-.2 

-.01 

-.2 

-.01 

50 

200 

90 

12. 

150 

-.1 

-.01 

-.1 

-.01 

50 

200 

90 

Notes: 

(1) 
(6) 

Depths  in  m 

Lengths  in  km 

Velocities  in  m/s,  +  to  north. 

Bottom  depth:  1  topography. 

Eddy  rotation:  CW  for  +  vel. 

Basin  rotation  angle  (y)  in  de^ 

-  to  south 
0  flat  bottom 
,  CCW  for  -  vel 
;rees  from  east 

.;  ne  = 

no  eddy 

similar  to   the  reference   state.     The  wave   structure  is  more   uniform  in   this  run, 
indicating  that  the  eddy  is  not  causative.  However,  it  does  aflect  the  wave  structure. 
c.   Experiment  No.  3  (Jet  Moved  Farther  Offshore) 

The  jet  and  the  eddy  were  moved  out  farther  away  from  the  coast  to  see  if 
wave  structures  were  being  generated  as  a  result  of  boundary  interactions  with  the 
eastern,  no-slip  boundary.  Figure  2.9  shows  the  output  fields  for  this  run.  In  general, 
the  results  are  very  comparable  to  the  reference  state.    It  looks  as  though  the  full 
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Figure  2.7     Experiment  #  I  (Hj  Increased  to  250  m) 
Model  Output  Height  Anomalies. 

Height  anomaly  fields  for  (a)  surface  at  dav  15  (cm),  (b)  surface  at  day  30  (cm),  (c) 
interface  at  day  15  (m),  and  (d)  interface  at  day  30  (m). 
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Figure  2.8     Experiment  #  2  (Eddy  Removed) 
Model  Output  Height  Anomalies. 

Height  anomaly  fields  for  (a)  surface  at  day  15  (cm|,  (b)  surface  at  day  30  (cm),  (c) 
interface  at  day  15  (m),  and  (d)  interface  at  day  30  (m). 
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growth  and  extent  of  the  waves  to  the  west  is  now  being  hindered  but  their  structure  is 
consistent.  Apparently,  the  boundaries  can  restrict  but  cannot  generate  the 
disturbances. 

d.  Experiment  No.  4  (  Topography  Removed) 

In  this  case,  topography  was  removed  and  a  flat  bottom  at  4500  m  depth 
utilized.  There  is  substantial  circumstantial  evidence  linking  filaments  to  coastal 
topographic  features  such  as  the  Mendocino  Escarpment.  By  removing  the 
topography,  it  is  hoped  that  it  can  be  determined  whether  bottom  variations  are 
causmg  the  filament  formation.  Figure  2.10  shows  the  output  fields  for  this  run.  The 
fields  are  nearly  identical  with  the  reference  run.  Apparently,  topography  is  not 
inducing  the  filament  formation  process  for  this  weak  lower  layer  flow. 

e.  Experiment  No.  5  (Increased  Lower  Layer  Velocity) 

In  this  run  the  lower  layer  velocity  was  increased  from  1  cm/s  to  5  cm/s 
flowing  southward,  which  is  considered  to  be  on  the  order  of  the  maximum  observed 
velocity  in  the  deep  water  south  of  Cape  Mendocino  (Stabeno  and  Smith,  1987).  This 
reduces  the  vertical  shear  but  increases  significantly  the  model  coupling  with  the 
bottom.  Although  the  previous  run  tended  to  discount  the  importance  of  the  bottom, 
it  is  important  to  see  what  possible  effects  it  might  have  in  a  worst  case  situation. 
Figure  2.11  shows  the  output  fields  over  the  domain.  The  results  are  nearly  identical 
with  the  reference  state  with  the  amplitudes  slightly  smaller.  This  indicates  that  the 
model  is  relatively  insensitive  to  southward  lower  layer  velocity. 

/.    Experiment  No.  6  (Lower  Layer  Flow  Direction  Reversed) 

In  this  experiment,  the  lower  layer  flow  was  reversed  so  that  it  was  flowing 
opposite  the  upper  layer  flow.  This  effectively  increases  the  vertical  velocity  shear  and 
examines  a  possible  configuration  of  the  CCS  where  there  is  a  deep,  weak  northward 
flow.  Figure  2.12  shows  the  output  fields  for  this  run.  The  wave  structure  is  similar 
to,  yet  noticeably  stronger  than,  the  reference.  This  would  indicate  the  wave  growth  to 
be  very  sensitive  to  the  vertical  velocity  shear. 

g.   Experiment  No.  7  (Zonal  Flow  Westward) 

In  order  to  perform  a  sort  of  check  on  the  model  dynamics,  the  basin  was 
not  rotated  into  its  CCS  position.  Since  most  instabiUty  studies  are  concerned  with 
zonal  rather  than  meridional  flows,  the  analytic  analysis  oi^  this  type  of  problem  is 
common  and  the  expected  dynamics  are  supported  by  observational  data.  With  y  =  0, 
the  jet  is  westward  and  the  eddy  rotates  as  before,  anticyclonically.    Figure  2.13  shows 
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Figure  2.9     Experiment  #  3  (Jet  Moved  Farther  Offshore) 
Model  Output  Height  Anomalies. 

Height  anomaly  fields  for  (a)  surface  at  day  15  (cm),  (b)  surface  at  day  30  (cm),  (c) 
interlace  at  day  15  (m),  and  (d)  interface  at  day  30  (m). 
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Figure  2.10    Experiment  #  4  (Topography  Removed) 
iModel  Output  Height  Anomalies. 

Height  anomaly  fields  for  (a)  surface  at  dav  15  (cm),  (b)  surface  at  day  30  (cm),  (c) 
interface  at  day  15  (m),  and  (d)  mterface  at  day  30  (m). 
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Figure  2.11     Experiment  #  5  (Increased  Lower  Layer  Velocity) 
iModel  Output  Height  Anomalies. 

Height  anomaly  fields  for  (a)  surface  at  dav  15  (cm),  (b)  surface  at  day  30  (cm),  (c) 
interface  at  day  15  (m),  and  (d)  interface  at  day  30  (m). 
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Figure  2.12    Experiment  #  6  (Lower  Layer  Flow  Direction  Reversed) 
Model  Output  Height  Anomalies. 

Height  anomaly  fields  for  (a)  surface  at  dav  15  (cm),  (b)  surface  at  day  30  (cm),  (c) 
interface  at  day  15  (m),  and  (d)  interface  at  day  30  (m). 
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the  output  from  this  simulation.  It  is  readily  apparent  that  there  are  some  type  of 
instabilities  in  the  flow  but  they  are  substantially  reduced  from  the  cases  where  there 
was  meridional  flow.  This  makes  sense  from  a  linear  instability  standpoint  and  will  be 
clarified  in  the  next  chapter. 

h.    Experiment  No.  8  {Zonal  Flow  Eastward) 

As  with  the  previous  experiment,  the  basin  was  not  rotated  into  its 
standard  CCS  orientation.  To  examine  the  model  response  for  a  jet  flowng  zonally 
eastward,  the  basin  was  rotated  180".  As  mentioned  before,  the  dynamics  of  this  type 
of  configuration  are  well  known  and  should  help  verify  the  model  responses.  Figure 
2.14  shows  the  output  fields  for  this  run.  Of  note  is  the  fact  that  the  flow  is  almost 
completely  stable  throughout  the  model  run.  Again,  this  makes  sense  from  a 
dynamical  standpoint  and  will  also  be  addressed  in  the  following  chapter. 
/.    Experiment  No.  9  {Eddy  and  Topography  Removed) 

In  experiment  #  2  the  eddy  was  removed  and  in  experiment  #  4 
topography  was  eliminated.  This  does  not  totally  resolve  their  involvement  with  the 
wave  structure  development.  A  requirement  for  baroclinic  instability  is  a  perturbation 
of  the  mean  flow.  It  is  conceivable  that  the  eddy  and  the  topography  could  be  causing 
the  jet  to  be  deflected  so  as  to  undergo  perturbation  growth  by  instability  theory.  To 
see  whether  this  is  true,  both  the  eddy  and  topography  were  removed  for  this  run.  As 
can  be  seen  in  Figure  2.15,  the  removal  of  both  the  eddy  and  topography  did  not 
hinder  at  all  the  development  of  the  wave  structures.  There  are  some  slight  differences 
from  the  reference  state  run,  but  they  are  so  small  that  it  would  be  difficult  to 
speculate  as  to  their  exact  connection  vn\h  the  altered  conditions  of  this  run  vice  the 
reference  one. 

j.   Experiment  No.  10  {Northwest  to  Southeast  Flow) 

As  can  be  seen  in  Figure  2.1,  the  coastline  south  of  Cape  Mendocino 
changes  orientation  from  north-south  to  northwest-southeast.  To  examine  the  effects 
of  this  orientation  on  the  flow  the  basin  was  rotated  120°  as  shown  in  the  figure  by  the 
dotted  box.  This  new  configuration  is  important  to  the  degree  of  applicability  the 
model  results  v^ll  have  to  the  CCS.  It  is  obvious  that  the  coastal  jet  cannot  continue 
straight  south  after  Cape  Mendocino  and  if  wave  structures  are  not  observed  with  the 
modeled  jet  going  southeast  it  is  doubtful  whether  the  meridional  flow  results  are  true 
representations  of  the  dynamics  of  the  CCS.  Figure  2.16  shows  the  model  response  to 
the  basin  rotation  of  120°.   The  eddy  was  removed  from  this  run  so  the  wave  structure 
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Figure  2.13     Experiment  #  7  {Zonal  Flow  Westward) 
Model  Output  Height  Anomalies. 

Height  anomaly  fields  for  (a)  surface  at  day  15  (cm),  (b)  surface  at  day  30  (cm),  (c) 
interlace  at  day  15  (m),  and  (d)  interface  at  day  30  (m). 
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Figure  2.14     Experiment  #  8  (Zonal  Flow  Eastward) 
Model  Output  Height  Anomalies. 

Height  anomaly  fields  for  (a)  surface  at  day  15  (cm),  (b)  surface  at  day  30  (cm),  (c) 
interface  at  day  15  (m),  and  (d)  interface  at  day  30  (m). 
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Figure  2.15    Experiment  #  9  (Eddy  and  Topography  Removed) 
Model  Output  Height  Anomalies. 

Heicht  anomaly  fields  for  (a)  surface  at  dav  15  (cm),  (b)  surface  at  day  30  (cm),  (c) 
interface  at  day  15  (m),  and  (d)  interface  at  day  30  (m). 
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changes  could  be  better  seen.  There  are  substantial  differences  from  the  reference  but 
they  are  clearly  in  terms  of  the  wave  structure  direction  and  wavelength  rather  than 
representing  a  totally  difTerent  type  of  phenomenon.  In  the  next  chapter,  the  reasons 
for  these  structural  differences  will  be  examined. 

k.    Experiment  No.  11  (V^-  Reduced  to  0.2  mjs) 

Up  until  now  all  of  the  model  runs  have  used  an  upper  layer  velocity  of  0.4 
m.  s.  This  may  be  representative  of  the  CCS  at  certain  times  but  by  no  means  is  the 
CCS  consistent.  As  with  the  previous  experiment,  appUcability  to  the  real  world 
demands  that  the  model  parameters  be  as  reaUstic  as  possible  in  both  value  and  range. 
The  velocity  is  reduced  to  0.2  m/s  thus  effectively  halving  the  vertical  shear.  Figure 
2.17  shows  the  results  of  this  alteration.  The  wave  structures  are  defmitely  present  in  a 
form  similar  to  the  reference  although  their  strength  is  considerably  less.  This  makes 
sense  from  an  instability  standpoint  as  the  shear  and  thus  available  potential  energy  of 
the  system  has  been  reduced.   This  \v\\\  be  discussed  in  the  next  chapter. 

/.    Experiment  No.  12  {V.-  Reduced  to  0.1  mjs) 

As  a  follow  on  to  the  previous  experiment,  the  upper  layer  velocity  is 
reduced  again  by  fifty  percent.  The  results  are  indicated  in  Figure  2.18.  There  is  some 
ver>-  slight  wave  development  but  it  is  considerably  less  than  either  the  previous  case  or 
the  reference  state.  This  would  intuitively  mean  that  the  reduction  of  the  upper  layer 
velocity  is  causing  the  system  to  approach  a  point  of  stability.  The  fact  that  the 
system  is  becoming  stable  with  a  shear  still  present  will  be  covered  also  in  the  next 
chapter. 
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Figure  2.16    Experiment  #10  (Northwest  to  Southeast  Flow) 
Model  Output  Height  Anomalies. 

Height  anomaly  fields  for  (a)  surface  at  dav  15  (cm),  (b)  surface  at  day  30  (cm),  (c) 
interface  at  day  15  (m),  and  (d)  interface  at  day  30  (m). 
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Figure  2.17    Experiment  #11  (Vj:  Reduced  to  0.2  m/s) 
iModel  Output  Height  Anomalies. 

Height  anomaly  fields  for  (a)  surface  at  dav  15  (cm),  (b)  surface  at  day  30  (cm),  (c) 
interface  at  day  15  (m),  and  (d)  interface  at  day  30  (m). 
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Figure  2.18     Experiment  #12  (Vj|  Reduced  to  0.1  m/s) 
Model  Output  Height  Anomalies. 

Height  anomaly  fields  for  (a)  surface  at  dav  15  (cm),  (b)  surface  at  day  30  (cm),  (c) 
interface  at  day  15  (m),  and  (d)  interface  at  day  30  (m). 
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III.  BAROCLINIC  ANALYSIS 

A.       SETTING  UP  THE  LINEAR  INSTABILITY  PROBLEM 

L  Instability  of  a  Non-Zonal  Flow 

Most  examinations  of  baroclinic  instability  start  with  the  assumption  of  a 
zonal  flow.  This  assumption  caimot  be  made  when  dealing  with  the  CCS  where  the 
basic  state  is  definitely  non-zonal  and  in  some  regions  completely  meridional.  To 
develop  an  analytic  model  for  the  CCS,  the  quasigeostrophic  equations  of  motion  will 
be  employed.  This  represents  a  divergence  from  the  model  dynamics  which  are 
operating  under  the  primitive  equations.  To  evaluate  whether  the  quasigeostrophic 
equations  are  reasonable  approximations,  the  Rossby  number,  Ro,  can  be  calculated. 
Using  U  =  0.4  m's,  f=0.94xio-^  s'^  and  L=Rj=18.4  km  yields  Ro  =  0.23  .  This 
indicates  that  nonlinearities  are  appreciable  but  the  quasigeostrophic  equations  siiould 
be  within  the  required  accuracy. 

A  non-zonal  flow  automatically  impUes  the  existence  of  an  external  forcing 
field  (Pedlosky,  1979).  In  this  case,  it  is  assumed  that  the  external  forcing  involved 
with  the  CCS  is  wind  stress  curl.  For  a  two  layer  system  without  bottom  friction  the 
governing  quasigeostrophic  potential  vorticity  equations  in  non-dimensional  form  irom 
Pedlosky  (1979)  are: 


=    [—^ ]     ;:    •    curl      t       ^^-^^ 


Notation  is  defined  in  the  Appendix. 

Following    the    derivation    by    Pedlosky   (1979)    for    the    phase    speed    of  a 
baroclinic  wave  in  a  non-zonal  flow,  consider  the  following  basic  state: 


Vj  =  Tj  =  VjX  -  u^y  (3.3) 
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M/,  =  'l\  =  0  --  (3.4) 

The  parameters  Uj  and  v^  arc  constant  velocities  in  the  upper  layer  in  tlic  x  and  y 
direction  respectively.  Note  that  the  basic  state  is  characterized  by  a  uniform  How 
which  does  not  vary  in  the  cross  stream  direction.  This  will  inherently  eliminate 
barotropic  instability  from  the  analytic  development  but  due  to  the  length  scales  of  the 
instabilities  observed  in  the  model  runs  (L>  >  Rj).  barotropic  instability  is  not  likely 
to  be  a  dominant  force  (Kamenkovich  et  ai,  19S6).  There  are  no  boundary  constraints 
on  the  basic  flow  ehminating  shear  induced  vorticity  which  is  present  in  the  numerical 
model.  Again,  it  is  felt  that  this  will  primarily  affect  barotropic  instability  modes.  The 
wall  effects  on  baroclinic  instability  are  coiisidered  negligible.  As  Mysak  ci  nl.  (19S1) 
pointed  out. 

.  .  .  inclusion  of  vertical  side  walls  at  the  edges  of  a  coastal  current  flowing  along 
a  continental  slope  does  not  significantly  affect  the  growth  rate,  frequency  or 
wavelength  of  the  maximum  unstable  wave  .... 

The  real  motivation  for  making  these  assumptions  is  that  the  analytic  modeling  of  a 
non-uniform,  non-zonal  flow  with  boundary  elTects  would  be  far  beyond  the  scope  of 
this  presentation  and  would  probably  not  yield  significantly  different  results.  Model 
runs  with  small  flows  in  the  lower  layer  in  both  a  north  and  south  direction  were  so 
similar  that  the  assumption  of  no  flow  in  the  lower  layer  should  be  consistent  with  the 
numerical  model. 

To   examine  the   stability  of  the   basic  state,  consider  perturbatior.s   of  the 
following  form: 

vj/,  =  T,  +  (p,(x,y.t)  (3.5) 


\\f,  =  T2  +  (P2(x.y,t)  (3.6) 

Substitution  into  3.1  and  3.2  and  linearization  viclds 


d  d  d  , 

I—  +   Uj—  +   Vj— ][  V-(pj-Fj((p,-(P2)I 
at  ^  dx  vy  1      I      I      ^ 


^(p,  ^(p,  (3.7) 


+    (P  +  FjU   )— •      +F  v,-^'=   0 
OK  '  c^y 
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[^][V'([^  -F,((p,-(p,)]    +    ((5-Fu   )V^-     -    F,v,^-     =   0     '      (?.8) 
c' t  i      ^      -  -^    '   c'x  -    '  c  y 

Note  that  the  wind  stress  curl  forcing  does  not  appear  in  the  perturbation  equations. 
As  Pedlosky  (1979)  describes  it, 

.  .  .  the  stability  of  basic  Hows  which  are  not  solutions  of  the  unforced  equations 
oi  motion  may  be  consistently  considered  in  the  context  o^  the  unforced 
perturbation  equations,  without  the  need  to  consider  explicitly  tlic  forces  required 
to  produce  a  basic  state. 

This  provides  the  basis  for  not  specifying  the  form  of  the  wind  stress  curl  in  3.1. 

Using  the  normal  mode  approach,  plane  wave  solutions  of  the  following  form 
are  assumed: 


v|/j^  =  Aj^  exp  i(kx  +  (y-(Tt)       n  =  1,  2  (3.9) 

Substitution  into  3.7  and  3.S  leads  to 

A,[  ((T-Ujk-v,e)(K2  +  F,  )    +    pk 

+    F,(u,k+v,C)]     -    A^F  (a-u  k-v  f)    =   0        ^^^'^'^^ 


A.laCK^+Fj)    +    pk   -     F2(Ujk+VjC)]     -    A^F2(T   =      0  (3.11) 

where  K"  =  k'  +  C^ 

A  coordinate  system  is  now  defined  as  shown  in  Figure  3.1.  The  angle  the 
perturbation  wave  vector  makes  with  the  x-axis  is  9  while  a  represents  the  angle  the 
shear  tlow  makes  with  the  x-axis. 
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k          u. 

Figure  3.1     Coordinate  System  Orientation. 

Based  on  the  flow  directions  of  the  CCS  within  the  numerical  model  domain  noted  in 
section  II. A,  a  can  take  on  only  one  of  two  values,  270°  or  300°.  The  coordinate 
system  allows  the  following  transformations  and  the  introduction  of  the  shear  velocity, 

Vs^ 


pk/K  =  P  cos  0 


(3.12) 


(u^k+Vj£)/K  =  Vgcosla-G) 


(3.13) 


Where  V^  =  (u^2  +  Vj2)l/2 


The  perturbation  wave  phase  speed  c  is  allowed  to  take  on  complex  values, 


c  =  (T/K  =  Cj.  ±  icj 


(3.14) 


After  dividing  3.10  and  3.11  by  K,  use  of  3.12,  3.13  and  3.14  results  in. 
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Ai[(c-V^cos(a-0))(K2+F^)    +      p^^^Q   ^    FjVgCosC  a-0)  ]       (3.15) 

-    A^FjC  c-V3Cos(u-0)  )    =   0 


AjcCK^+F^)    +    PcosB    -    F^VgCosCa-O)]     -   AjF^c    =   0         (3.16) 

Nontrivial  solutions  for  A^  and  A2  exist  only  if  the  deterniinant  of  their  coefficients  in 
3.15  and  3.16  is  equal  to  zero.  Since  the  basic  intent  of  this  analytic  development  is 
specifically  to  support  analysis  of  the  model  output  and  the  CCS  it  is  more  practical  at 
this  point  to  switch  to  dimensional  form.  Using  the  notation  of  Kamenkovich  ei  al. 
(1986),  the  zero  determinant  condition  leads  to  the  phase  speed  of  the  perturbation 
wave  (Equation  3.17). 

V3K^(K-+2AOcos(a-0)    -     P(  2K-+A,+AOcosO 

c    = ! :: 

2K^(K^+Ai+A,) 

[  (Aj  +  A2)-p^cos^e   +   2pV3K-*(Aj-A2)cosecos(a-e)  (3.17) 

-    K'*V^'(4A,A,-K'^)cos-(a-e)]^/2 

±    2 i—:^ 

2K-(K-+A,+A2) 

Where  Aj  =  Fj/L^,  A2=F2/L^,  P  =  Pq  and  all  velocities  from  before  multiplied  by 
scaling  velocity  U. 

2.  Critical  Velocity  and  Growth  Rate 

When  the  radicand  in  3.17  is  equal  to  zero,  the  wave  is  at  the  transition  point 
between  stability  and  instability.  At  this  point,  V^  can  be  solved  for.  This  results  in  an 
expression  for  the  critical  shear  velocity  for  instability,  V^^. 


PK2(A,-A,)cose    ±    2P(A,A,)V2[(A.+A,)--K'*]^/2^Q3Q 

V        -  ' 


sc  K2(4AjA2-K'*)cos(a-e) 


(3.18) 
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Utilizing  tiie  two  basic  state  flow  directions  to  match  the  CCS  (a  =  270°, 
300°)  and  using  model  values  of  A^  =  2.94  x  10"^  m'"^,  A2=  1.02  x  10'^*^  m'^  and 
P=2.0x  10  m  s'^  leaves  V^^  a  function  of  only  K  and  G.  Figures  3.2  and  3.3  show 
the  marginal  stabiUty  curves  for  a  of  270°  and  300°  respectively  for  various 
perturbation  wave  directions,  0.  Note  that  curves  are  not  shown  for  cases  where  the 
wave  direction  is  perpendicular  to  the  flow  direction  since  by  3.18  this  results  in  infinite 
shear  velocities  and  no  instability  is  possible.  The  form  of  the  marginal  stability  is 
determined  by  the  quadrant  that  G  is  in.  The  unstable  region  is  shaded  and  represents 
the  area  where  the  shear  velocity  is  greater  than  required  for  instability.  The  excess 
shear  then  becomes  the  energy  source  to  provide  a  positive  growth  rate  in  those 
regions.  The  full  implications  of  a  meridional  flow  and  a  meridional  wave  direction  are 
indicated  by  Figures  3.2b  and  3.2e..  The  required  shear  flow  goes  to  zero  and  thus  any 
non-zero  shear  will  result  in  instability  as  long  as  the  wavenumber  is  to  the  left  of  the 
vertical  asymptote  corresponding  to  the  condition, 

K  <  (^A^A^)^/"^  (3.19) 

which  comes  directly  from  Equation  3.18. 

With  the  shear  conditions  for  instability  defined,  the  next  step  is  to  find  the 
magnitude  of  the  growth  rate  in  the  unstable  regions.  Growth  rate  values  will  provide 
the  necessary  information  to  determine  the  characteristics  of  the  fastest  growing  wave. 
For  an  unstable  wave  to  grow,  the  radicand  of  3.17  must  be  negative  and  the  growth 
rate,  Kc-  becomes, 

[KVs2(4AjA2-K*)cos2(a-e)- 
_  -   2pVgK^(A^-A2)cosecos(a-9)    -   ( Aj-^AQ^p-cos-Ql  V2 

2K(K^-^A,*A,)  (320) 


As  before,  model  values  for  A^,  A2  and  p  along  with  flow  angle  a  reduces  the  growth 
rate  dependence  to  just  V^,  9  and  K.  Figures  3.4  and  3.5  show  the  growth  rates  for  a 
=  270°  and  300°  respectively  for  velocity  shears  of  0.4,  0.3,  0.2  and  O.l  m/s.  As 
expected,  the  growth  rate  magnitude  is  a  strong  function  of  the  shear  velocity  since 
that  provides  the  energy  input  into  the  growing  wave.    The  maximum  growth  rate, 
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Figure  3.2    Marginal  Stability  Curves  for  0  =  270". 
Curves  for  9  of  (a)  45^  (b)  90^  (c)  I35^  (d)  225^  (e)  270"  and  (0  315' 
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Figure  3.3     Marginal  Stability  Curves  for  a  =  300°. 
Curves  for  0  of  (a)  75°,  (b)  120  **  (c)  165^  (d)  255^  (e)  300°  and  (0  345' 
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Figure  3.4    Growth  Rates  for  a  =  270°. 
Growth  rates  (  x  10"^  s"^)  for  velocity  shears  of  (a)  0.4,  (b)  0.3,  (c)  0.2  and  (d)  0.1  m/s. 
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Figure  3.5    Growth  Rates  for  a=  300  . 
Growth  rates  (  x  10"^  s"^)  for  velocity  shears  of  (a)  0.4,  (b)  0.3,  (c)  0.2  and  (d)  O.I  m/s. 


49 


which  would  correspond  to  the  fastest  growing,  wave,  occurs  at  difTerent  K-9 
intersections  in  each  difTerent  shear  flow  case.  This  means  the  shear  flow  not  only 
affects  the  magnitude  of  the  growth  rate  but  also  the  direction  and  wavenumber  of  the 
perturbation  wave.  As  a  result,  any  variation  in  the  external  forcing  which  is 
producing  the  shear  flow  would  cause  major  alterations  to  the  fastest  growing  wave. 
Another  interesting  point  is  that  the  growth  rates  for  0  =  270*  are  higher  in  all  shear 
conditions  than  for  0=300*  although  the  G-K  plane  positive  growth  rate  areal  extent 
is  smaller.  From  the  marginal  stability  curves,  it  is  apparent  that  a  meridional  flow 
could  have  perturbations  for  which  any  non-zero  shear  would  support  instability.  As  a 
result,  the  additional  claim  that  meridional  flow  yields  the  largest  growth  rate  for  a 
given  shear  can  be  made. 

B.       APPLICATION  OF  ANALYTIC  RESULTS  TO  MODEL  OUTPUT 

The  instability  problem  has  now  been  formulated  to  approximate  the  numerically 
modeled  situation.  The  next  step  is  to  examine  the  instabilities  present  in  the  model 
output  fields  as  noted  in  section  n.B.2  in  terms  of  the  previous  section's  results.  If 
there  are  no  inconsistencies,  then  the  instabilities  are  most  likely  due  to  baroclinic 
instability  of  the  basic  fiow. 

From  section  II. B. 2,  it  is  obvious  that  there  is  a  similarity  of  the  output  fields 
with  respect  to  apparent  wave-like  structures.  Not  only  do  the  same  features  appear  in 
most  of  the  runs,  but  their  shape,  direction  and  magnitudes  are  also  similar.  Table  3 
shows  the  average  values  of  relevant  parameters  from  the  model  runs  (where  observed) 
for  days  10,  20  and  30  for  flow  direction  of  a=270°.  Table  4  shows  the  same 
parameters  for  flow  direction  of  a=  300*.  The  wave  data  for  both  tables  was  taken 
from  model  runs  where  the  shear  flow  was  close  to  0.4  m/s  ( ±  0.01  m/s). 


TABLE  3 

WAVE  PARAMETERS  FROM  MODEL  OUTPUT  FIELDS  (a=270°) 

DAY 

G  (deg. 

)                K{m-b 

L(km) 

10 

215 

2.0  X   10"^ 

310.0 

20 

249 

2.6  X   10'^ 

240.0 

30 

255 

3.1   X   10"^ 

200.0 

50 


TABLE  4 

WAVE  PARAMETERS  FROM  MODEL  OUTPUT  FIELDS  (a  =  300") 

DAY 

6  (deg.) 

K(m-b 

L(km) 

10 

238 

2.3  X   10-^ 

275.0 

20 

260 

3.0  X   10-^ 

209.0 

30 

269 

3.2  X   10"^ 

198.0 

These  days  were  selected  because  they  represent  the  appearance  (day  10),  the 
growth  (day  20)  and  the  maturity  (day  30)  of  the  waves.  By  mapping  the  wave  data 
from  Tables  3  and  4  onto  instability  plots  from  section  in.A.2,  an  evaluation  of  their 
relationship  to  the  analytically  derived  baroclinic  wave  characteristics  can  be  made. 

Marginal  stability  curves  are  shown  in  Figure  3.6  for  the  flow  angles  and  wave 
angles  noted  in  Tables  3  and  4.  The  wavenumbers  for  those  angle  values  are  indicated 
and  marked  by  the  day  the  data  is  associated  with  at  its  intersection  with  the  ¥^  =  0.4 
m/s  line.  Of  note  is  the  fact  that  all  of  the  plots  show  the  waves  within  the  unstable 
region  for  baroclinic  instability.  Additionally,  as  the  waves  shift  in  wavenumber 
towards  shorter  wavelengths,  their  shear  excess  increases  to  its  maximum  where  the 
wave  meets  the  asymptote  corresponding  to  Equation  3.19. 

Growth  rate  plots  are  shown  in  Figure  3.7  for  velocity  shears  of  0.4  m/s.  The 
parameter  positions  corresponding  to  days  10,  20  and  30  from  Tables  3  and  4  are 
indicated  to  show  the  path  the  waves  took  as  they  developed.  From  Figure  3.7  it  is 
apparent  that  the  model  output  waves  fall  in  the  regions  of  positive  growth  rate  but  do 
not  fall  on  or  near  the  peak  growth  rate  which  would  correspond  theoretically  to  the 
position  of  the  fastest  growing  wave.  Also,  the  growth  rate  is  low  at  day  10,  reaches  a 
maximum  at  day  20  and  is  decreasing  by  day  30.  It  must  be  kept  in  mind  that  these 
growth  rates  represent  the  analytic  values  given  the  wavelength  and  wave  direction  of 
the  model  output  waves.  They  are  not  necessarily  the  same  as  the  growth  rates  of  the 
actual  waves  in  the  model  output.  The  basic  effect  of  the  growth  rate  is  to  increase  the 
amplitude  of  the  growing  wave  in  the  following  manner: 

A  =  Age^M  (3.21) 
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Figure  3.6    Marginal  Stability  Curves  for  Model  Output  Wave  Characteristics. 

Curves  for  a  =  270*  anil  0=  (a)  215",  (b)  249°,  (c)  255"  and  for  0=300"  and  9=  (d) 
238\(e)260   ,(f)269    . 
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Figure  3.7    Growth  Rate  Positions  of  Model  Output  Waves  for  Shear  Velocity  of  0.4  m/s. 

Parameter  data  from  (a)  Table  3  and  (b)  Table  4  ploted  on  growth  rate  (  x  10*"  s"^) 
fields  from  Figures  3. 4. a  and  3. 5. a. 

where  A  is  the  time  varying  amplitude  and  Ag  is  the  initial  amplitude.    Rearranging 
Equation  3.21, 


£n(A  -  An)  =  kcjt 


(3.22) 


yields  an  expression  of  the  form,  y=  mx  where  m  represents  the  slope  (kcj).  By  plotting 
the  natural  log  of  the  amplitude  difference  from  the  initial  versus  time,  the  growth  rate 
can  then  be  determined  by  the  slope  of  the  curve.  Figure  3.8  shows  the  amplitude 
changes  of  the  reference  state  waves  during  the  model  run. 

Immediately  apparent  is  the  fact  that  the  growth  rate  varies  with  time.    It  is  low 
at  first  and  grows  to  a  maximum  and  then  decreases.   The  information  from  Figure  3.8 
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Figure  3.8     Reference  State  Wave  Amplitude  Growth. 

can  be  easily  converted  into  a  plot  of  growth  rate  versus  time.  This  is  shown  in  Figure 
3.9.  The  temporal  variation  of  the  growth  rate  is  more  clear  and  shows  similar  form  to 
the  growth  rate  progression  observed  in  Figure  3. 7. a. 

Experiments  7  and  8  involved  zonal  flows.  These  were  not  meant  to  simulate  the 
CCS  but  to  evaluate  the  model  response  against  theor>'.  For  a  zonal  flow,  the  wave 
growth  must  offset  one  hundred  percent  of  the  stabilizing  planetary  vorticity  gradient, 
whereas  for  a  meridional  flow  the  gradient's  effects  are  zero  due  to  the  PcosG  term. 
Figure  3.10  shows  the  marginal  stability  curves  for  an  eastward  and  westward  zonal  jet. 
The  curves  are  not  for  any  specific  0  because  the  0  dependence  cancels  out  of  equation 
3.18  for  a  zonal  flow.  Immediately  apparent  is  the  fact  that  the  easterly  jet  should  be 
more  stable  than  the  westerly  jet  since  it  has  a  much  higher  minimum  shear  velocity. 
This  agrees  with  the  experimental  results  where  the  easterly  jet  remained  stable  through 
the  entire  run  and  the  westerly  jet  developed  slight  instabilities.    Based  on  the  curves 


54 


Figure  3.9    Growth  Rate  Variations  of  Reference  State  Waves. 

for  the  zonal  flow  coupled  with  the  fact  that  the  shear  velocity  was  0.4  m/s,  one  might 
expect  instability  wave  growth  for  both  zonal  flows  similar  to  the  reference  run 
although  this  was  not  observed.  The  increase  in  the  required  V^^  for  the  zonal  flows 
due  to  the  P  effect  cuts  down  on  the  excess  shear  available  for  growth  and  apparently 
this  excess  is  not  quite  great  enough  for  wave  growth.  This  points  to  frictional  eflccts 
and  will  be  discussed  in  the  next  section. 


C.       FURTHER  ASPECTS  OF  THE  INSTABILITY  PROBLEM 

The  assumptions  made  to  support  the  application  of  baroclinic  instability  theory 
to  the  numerical  model  results  were  necessary  for  any  quantitative  analysis  to  be  made. 
Generally,  there  was  very  good  agreement  between  linear  instability  theory  and  the 
model  output  which  supports  the  contention  that  baroclinic  instability  is  in  fact  the 
dominant  dynamical  process  in  the  model  and  in  the  CCS.  Nevertheless,  there  were 
discrepancies  between  the  model  simulations  and  the  analytic  predictions  and  these 
need  to  be  addressed, 
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Figure  3.10    Marginal  Stability  Curves  for  Zonal  Flow. 
Curves  are  for  a=  270°  and  flows  to  (a)  east  and  (b)  west. 

1.  Effects  of  Friction 

Without  friction,  the  vertical  phase  shift  of  a  baroclinic  wave  opposite  the 
vertical  velocity  shear  is  an  indication  that  c->0  and  the  wave  should  grow  (Pedlosky, 
1979).  However,  when  frictional  forces  are  introduced,  the  wave  can  take  on  the  same 
phase  shift  opposite  the  velocity  shear  and  still  be  neutral.  Instead  of  supplying  energy 
for  wave  growth,  the  extracted  energy  is  spent  offsetting  dissipation.  Thus,  in  order  to 
grow  in  the  presence  of  friction,  the  baroclinic  wave  must  have  a  greater  phase  shift 
than  would  be  the  case  without  friction.  Since  there  is  a  maximum  phase  shift 
associated  with  maximum  energy  extraction,  the  growth  rate  of  a  baroclinic  wave  with 
friction  should  be  less  than  that  for  a  wave  without  friction. 

In  Figures  3.4  and  3.5  it  was  observed  that  changes  in  the  shear  velocity 
resulted  in  maximum  growth  rate  value  and  0-K  positional  changes.  Frictional  forces 
would  tend  to  reduce  the  available  energy  for  growth  which  is  equivalent  to  a  reduction 
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in  the  shear  velocity.  For  this  reason,  the  numerical  output  waves  could  not  be 
expected  to  fall  on  the  maximum  growth  rate  parametric  positions  predicted  by  the 
non-frictional  analytic  development.  Another  complicating  factor  is  the  way  in  which 
friction  was  parameterized  in  the  numerical  model.  By  using  a  A^V^V  form,  friction  is 
a  strong  function  of  the  velocity  field.  Since  the  velocity  field  in  the  domain  changed 
drastically  as  the  baroclinic  waves  developed,  the  frictional  forcing  must  also  have 
changed  drastically.  Thus,  the  frictional  alteration  of  the  maximum  growth  rate  9-K 
position  would  take  on  temporal  variability.  This  might  explain  one  reason  the  model 
output  waves  were  observed  to  vary  in  wavenumber  and  direction  through  the  model 
runs. 

2.  Non-linear  Instability  Effects 

In  linear  instability  theory  the  perturbation  wave  is  unbounded  as  long  as  it 
meets  the  necessar>'  conditions  for  instability  which  are  defined  in  terms  of  the  basic 
state.  Obviously,  this  is  not  realistic  as  the  waves  grow  to  finite  amplitude.  Non-linear 
theory  explains  how  advection  affects  the  wave  growth. 

As  Pedlosky  (1979)  pointed  out,  it  is  the  basic  state  in  the  absence  of  the 
perturbation  wave  that  must  be  used  as  criterion  for  instability  vice  any  instantaneous 
fiow  profile.  Variations  in  the  flow  profile  as  the  wave  develops,  however,  does  have 
an  elTect  on  the  wave.  When  the  wave  reaches  a  point  where  the  instantaneous  mean 
fiow  would  not  support  instability,  the  wave  growth  rate  is  at  its  maximum  and  will 
start  to  decline.  Although  declining,  the  growth  rate  is  still  positive  and  the  wave  is 
still  growing  even  though  the  instantaneous  mean  fiow  might  support  no 
instability. The  change  in  growth  rate  is  accompanied  by  a  reduction  of  the  phase  shift 
of  the  wave  in  the  vertical  until  growth  is  stopped  as  the  phase  shift  reaches  zero.  The 
process  does  not  stop  here  but  continues  into  the  phase  orientation  that  causes  wave 
decline.  As  with  the  growth,  the  decline  will  hit  a  maximum  and  reverse  and  the  whole 
process  of  growth  and  decay  will  oscillate  back  and  forth  (Pedlosky,  1979).  Thus,  the 
variations  in  the  wave  growth  rate  observed  in  Figure  3.9  and  inferred  from  Figure  3.7 
could  well  be  indicative  of  non-linear  effects  present  in  the  model  and  presumably,  the 
CCS. 

When  considering  the  effects  of  friction,  variations  in  the  available  energy  to 
the  wave  were  correlated  to  variations  in  the  parametric  characteristics  of  the 
maximum  growth  wave.  Non-linear  theory  predicts  similar  variations  in  the  available 
energy  due  to  the  changes  in  the  mean  fiow  brought  about  by  the  perturbation  wave. 
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Hence,  it  is  another  possible  causative  factor  for  the.  model  output,  since  the  model  is 
non-linear,  to  exhibit  variability  in  baroclinic  wave  wavenumber  and  direction  over 
time. 

3.  Mixed  Instability  Effects 

Most  treatments  of  the  instability  problem  are  either  focused  on  baroclinic  or 
barotropic  instability.  In  the  CCS  and  the  numerical  model,  flows  have  both 
horizontal  and  vertical  shears.  It  was  assumed  in  the  analytic  development  of  the 
baroclinic  instability  problem  in  section  III.A.l  that  barotropic  instability  was  not  the 
dominant  instability  process  based  on  the  scales  of  the  observed  disturbances. 
However,  it  may  be  that  barotropic  instability  does  have  a  slight  influence  on  the 
baroclinically  unstable  waves. 

Pedlosky  (1979)  showed  that  for  a  system  with  vertical  and  horizontal  velocity 
shears,  the  characteristics  of  the  most  unstable  wave  can  only  be  derived  with 
exhaustive  calculation  given  a  detailed  spatial  distribution  of  velocity  and  potential 
vorticity.  Changing  the  parameters  of  the  flow  can  have  a  strong  inpact  on  the  type  of 
instability  present  and  the  wavelength.  It  is  believed  that  throughout  the  numerical 
simulations,  baroclinic  instability  remained  the  dominant  instability  process  but 
barotropic  instability  may  have  altered  the  characteristics  of  the  fastest  growing  wave. 

4.  Baroclinic  Instability  with  Complex  Wavenumbers 

In  developing  the  baroclinic  instability  problem,  a  particular  perturbation 
form  was  assumed.  This  form  is  common  to  nearly  all  studies  of  this  type  of  instability 
and  is  shown  by  Equations  3.9  and  3.14.  The  crucial  point  is  the  selection  of  a  real 
wavenumber  and  a  complex  phase  speed. 

Hogg  (1976)  examined  analytically  and  numerically  the  baroclinic  instability 
process  in  a  zonal  flow  where  the  perturbation  form  allowed  both  complex  phase  speed 
and  complex  wavenumber.  In  general,  his  results  indicated  a  wavenumber  cutoff 
between  spatially  growing  (k->0)  and  temporally  growing  (c->0)  waves.  Temporally 
growing  waves  were  confined  to  long  wavelengths  while  spatially  growing  waves  were 
confined  to  the  shorter  wavelengths. 

In  a  study  of  Gulf  Stream  meanders,  Watts  and  Johns  (1982)  used  statistical 
analysis  of  cross  stream  data  and  derived  values  for  both  temporal  and  spatial  growth 
rates  of  the  meanders.  Contrary  to  Hogg  (1976),  Watts  and  Johns'  (1982)  results 
indicated  that  spatial  and  temporal  growth  rates  can  occur  in  similar  wavelength 
regimes.  Figure  3.11  shows  the  growth  rate  versus  wavenumber  relation  for  both 
spatial  and  temporal  growth  found  by  Watts  and  Johns  (1982). 

58 


' 

a* 

GROWTH 

- 

^-•-^ 

RATE 

' 

(I0'2  km-')  04 

■ 

1 

1                     1 

02 

GROWTH 
RATE 
(d-')        °' 


WAVENUMBER     k,    (IQ-^km'') 


Figure  3. 1 1     Spatial  and  Temporal  Wave  Growth  of  Gulf  Stream  Meanders 

(Watts  and  Johns,  1982). 

There  is  evidence  both  analytically  and  observationally  that  spatial  growth 
may  be  involved  in  the  baroclinically  unstable  wave.  Although  Hogg  (1976) 
concentrated  on  spatial  growth  (kj>0)  and  Watts  and  Johns  (1982)  observed  only 
temporal  and  spatial  growth  (C|>0,  k->0),  there  is  nothing  in  the  theory  that  prevents 
a  combination  of  temporal  growth  and  spatial  decline  (Cj>0,  k.j<0).  This  might  also 
explain  the  variation  in  the  observed  wavelengths  in  the  model  output  from  long  to 
short  while  the  waves  were  growing  in  time. 
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IV.  DISCUSSION  AND  CONCLUSIONS 

A.       ANALYSIS  OF  RESULTS 

I.  Comparison  of  Results  with  Observational  Data 

As  mentioned  in  the  introduction,  filament  structure,  spatial  distributions  and 
temporal  variabilities  are  well  known.  In  order  to  validate  the  results  of  this  study, 
there  should  be  general  agreement  with  the  existing  information  and  if  discrepancies 
are  present  there  should  be  some  theoretical  basis  for  the  differences.  In  Table  5 
values  are  shown  for  the  range  and  magnitude  of  the  observed  versus  the  modeled 
characteristics.  The  value  for  the  observed  wavelength  and  direction  is  from  Ikeda  and 
Emery  (1984)  and  the  observed  offshore  extent  is  from  Brink  and  Hartwig  (1985).  The 
rest  of  the  observed  values  are  from  Flament  et  al.  (1985). 
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TABLE  5 

FILAMENT  CHAR.\CTERISTICS,  OBSERVED  VS.  MODELED 

PARAMETER 

OBSERVED 

MODELED 

Wavelength 

100-400  km 

198-310  km 

Wave  Direction 

-270° 

215-269° 

Off-shore  Extent 

100-500  km 

100-800  km 

Maximum  Off-shore  Velocity 

-0.5  m/s 

0.5  m/s 

Maximum  On-shore  Velocity 

-0.35  m/s 

0.4  m/s 

Growth  rate  (e-folding  lime) 

— 1-2  weeks 

-lOd 

It  would  seem  that  the  only  parameter  that  deviated  significantly  was  the  wave 
direction.    The  reason  for  this  will  be  discussed  in  the  next  section.    Otherwise  the 
agreement  strongly  supports  the  contention  that  the  modeled  phenomena  are  in  fact 
representative  of  filaments  in  the  CCS. 
2.  The  Filament  Formation  Process 

The    filament    formation    process    is    quite    complex.     Although    baroclinic 
instability  is  clearly  the  dominant  mechanism,  there  are  other  forces  at  work  that 


cannot  be  neglected. 
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In  the  first  phases  of  growth,  the  wave  is  .affected  by  the  east-west  horizontal 
velocity  shear.  Close  to  the  jet,  the  wave  front  is  advected  strongly  to  the  south. 
Farther  offshore,  the  jet  effect  on  the  wave  propagation  decreases  due  to  the  mean 
velocity  gradient.  This  causes  the  wave  vector  direction  to  be  deflected  from  the 
seemingly  more  favorable  southerly  direction  to  a  near  westerly  direction.  As  the  wave 
gains  energy  from  the  vertical  shear  of  the  mean  flow,  advection  of  the  mean  flow  in  a 
zonal  direction  starts  to  break  down  the  east-west  horizontal  velocity  shear  and  the 
wave  direction  rotates  towards  its  more  optimum  direction  and  its  wavelength  shortens 
to  maximize  the  growth  rate.  Non-hnearities  become  dominant  as  the  wave  reaches 
such  a  magnitude  that  the  instantaneous  mean  flow  is  stable  and  the  wave  growth 
starts  to  decline.  This  is  as  fir  as  the  waves  were  modeled,  but  it  would  be  expected 
that  the  waves  would  go  into  an  oscillatory  growth-decline  mode  until  damped  by 
friction. 

Obviously,  this  whole  analysis  is  contingent  on  a  constant  wind  stress  forcing 
for  periods  on  the  order  of  months  which  is  hardly  realistic.  It  can  be  expected  that 
constant  forcing  may  persist  for  weeks  at  a  time  in  which  case  the  initial  growth  and 
development  holds  true. 

From  an  observational  standpoint,  the  filaments  observed  are  nearly  always 
oriented  perpendicular  to  the  flow  indicating  a  southerly  wave  direction.  This  does  not 
mean  that  the  described  theory  of  development  is  necessarily  flawed,  it  just  means  that 
by  the  time  the  waves  are  observable,  they  have  probably  already  twisted  around  into 
their  more  optimum  configuration. 

B.       IMPORTANCE  OF  RESULTS 

There  have  been  numerous  studies  of  zonal  instability  processes  and  a  few  of 
meridional  situations  but  there  is  a  large  and  distinct  scarcity  of  studies  of  flows  in- 
between.  When  one  considers  currents  such  as  the  Gulf  Stream,  the  Kuroshio  and,  of 
course,  the  CCS  which  have  flows  that  are  not  always  east-west  or  north-south,  it 
becomes  evident  that  more  attention  should  be  placed  on  this  discrepancy. 

As  demonstrated,  the  behavior  of  a  flow  that  is  non-zonal  is  dramatically 
different  from  a  zonal  flow.  In  a  zonal  flow,  the  energy  releasing  horizontal 
temperature  gradient  and  the  stabilizing  planetary  vorticity  gradient  are  always  working 
against  each  other.  In  a  meridional  flow,  energy  release  can  take  place  without  any 
component  of  P  acting  on  it;  hence,  any  shear  above  the  dissipation  level  may  produce 
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instability.  For  the  in-between  case,  the  wave  must  assume  a  balance  between  the 
maximization  of  energy  release  and  the  dampening  effects  of  the  planetary  vorticity 
gradient. 

C.  RECOMMENDATIONS  FOR  FUTURE  STUDIES 

There  are  several  aspects  of  this  modeling  work  that  could  be  modified  to  better 
simulate  the  CCS.  A  better  resolution  model  would  help  significantly  to  bring  the 
barotropic  instability  process  into  view.  This  is  particularly  important  to  the  filament 
situation  since  observational  evidence  has  shown  instabilities  along  the  filament  frontal 
region  of  the  scale  of  R^^  that  are  claimed  to  be  the  result  of  barotropic  instability 
(Flament  ei  ai,  1985).  A  model  with  more  layers  would  give  the  vertical  resolution 
necessan.-  to  properly  simulate  the  CCS  in  all  of  its  configurations,  particularly  its  late 
summer  profile  of  a  strong  southward  surface  jet  and  a  poleward  undercurrent  all 
above  a  relatively  quiescent  deep  water  region.  The  addition  of  some  temporal 
variability  to  the  coastal  jet  strength  would  help  to  better  understand  wind  stress 
variation  elTects  on  existing  filaments  and  the  growth  of  future  filaments.  Continued 
intensive  in  situ  and  satellite  gathering  of  oceanographic  data  is  vitally  important  to 
understanding  the  physical  processes  necessary  to  better  reproduce  the  observed 
filament  structures  numerically.  To  thoroughly  determine  if  baroclinic  instability  is 
occurring  in  the  CCS.  an  elaborate  experiment  would  have  to  be  designed  to  resolve 
spatially  and  temporally  the  mesoscale  motions.  An  experiment  comparable  to  the 
POLYMODE  experiment  in  the  North  Atlantic  would  be  necessary. 

D.  CONCLUSIONS 

Baroclinic  instability  is  a  prime  generating  force  for  filaments  formed  south  of 
Cape  Mendocino  in  the  CCS.  Barotropic  instability  along  with  non-linearities, 
frictional  effects  and  possible  spatial  wave  growth  may  help  to  modify  the 
characteristics  of  the  growing  baroclinically  unstable  waves. 

Filament  growth  is  strongly  dependent  on  the  vertical  velocity  shear,  the  vertical 
structure  and  the  mean  flow  direction.  Any  attempt  to  model  filaments  must  carefully 
parameterize  these  quantities  if  accurate  and  thus  useable  results  are  to  be  obtained. 
Indeed,  for  any  numerical  simulation  of  currents  with  unknown  instabilities,  their 
values  are  crucial  to  correct  simulations.  Improper  understanding  of  their  spatial  and 
temporal  characteristics  can  seriously  affect  the  validity  of  the  model  results. 
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APPENDIX 
SYMBOLS  AND  NOTATION 

A^  Laplacian  lateral  friction  coefficient 

A^  Amplitude  coefficient  for  normal  mode  solution  for  n     layer 

c  Phase  speed  of  perturbation  wave 

D(x,y)  Variable  depth  of  topography 

Dj  Upper  layer  thickness 

F-  Froude  number  for  i'-"  layer  =  fg^L^'  g'H^ 

fg  Coriolis  parameter  for  mean  latitude 

g  Gravitational  acceleration 

g'  Reduced  gravitational  acceleration  =  giPj-fiiJ/Pi 

h:  Instantaneous  laver  thickness 
1 

H-  Upper  (i=  1)  and  lower  (i=  2)  layer  mean  thickness 

k  Unit  vector  in  the  z  direction 

k  Perturbation  wave  number  in  the  x  direction 

t  Perturbation  wave  number  in  the  y  direction 

K  Total  wave  number  for  perturbation  wave  =  k^  +  £^ 

L  Horizontal  scale  length 

p  I  Pressure  in  the  upper  layer  =  g(hj  +  h2  +  d) 

p2  Pressure  in  the  lower  layer  =  p^-g'h| 

R^  First  internal  Rossby  radius  of  deformation  =[g'H^H2/fQ  (Hj +  H2)] 

Rq  Rossby  number  =  U/fL 

U  Scale  velocity 

u-,v-  Velocities  in  the  x  and  y  directions 

U-,V-  Transport  in  the  x  and  y  directions 

V^  Shear  velocity  =(u,^  + v^^)*''^ 

Vg^  Critical  shear  velocity  for  instability 

x,y  Cartesian  coordinates  directed  N  and  W  respectively 

a  Angle  mean  flow  vector  makes  with  x  axis 

Pq  Variation  of  Coriolis  parameter  with  latitude  =  df^'dy 

P  Scaled  variation  of  Coriolis  parameter  with  latitude        =  PqL^/U 

Y  Basin  rotation  angle  from  east 

Ax,Ay  Grid  spatial  resolution 
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At  Model  time  increment 

9  Angle  perturbation  wave  vector  makes  with  x  axis 

A-  Dimensionalized  Froude  number  for  i^^  layer  =  F:'  L^ 

p-  Density  in  i^^  layer 

Pg  Mean  density 

(7  Frequency  of  perturbation  wave 

T  Wind  stress  vector 

(p-  Perturbation  stream  function  for  i      layer 

T-  Basic  state  stream  function  for  i     layer 

V|/-  Variable  stream  function  for  i^  layer 

V  Gradient  operator  =d;dx  +  d.'dy 

V^  Laplacian  operator  =d^idy?  +  c^jdy'^ 
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